Draft version February 19, 2013 

Preprint typeset using K^T^^X style emulateapj v. 5/2/11 



THERMODYNAMICS OF THE COMA CLUSTER OUTSKIRTS 
A. SiMiONESCU^'^, N. Werner^'^, O. Urban^'^'*, S. W. Allen^-^'*, A. C. Fabian^, A. Mantz"'-", K. Matsushita'', 

P. E. J. NULSEN*, J. S. SANDERS^■^ T. SASAKI^ T. SATO^ Y. TAKEi"', S. A. WALKER-* 

iRIPAC, Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA 
^Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060, USA 
3SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 
^Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK 
^Kavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 
^Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637-1433, USA 
^Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan 
^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 
®Max Planck Institute for Extraterrestrial Physics, Giessenbachstr 1, 85748 Garching, Germany and 
^"Institute of Space and Astronautical Science (ISAS), JAXA, 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawra, 252-5210 Japan 

Draft version February 19, 2013 



ABSTRACT 

We present results from a large mosaic of Suzaku observations of the Coma Cluster, the nearest 
and X-ray brightest hot (~8 keV), dynamically active, non-cool core system, focusing on the thermo- 
dynamic properties of the ICM on large scales. At azimuths not aligned with an infalling subcluster 
towards the southwest, our measured temperature and X-ray brightness profiles are in broad agree- 
ment, with the temperature decreasing from about 8.5 keV at the cluster center to about 2 keV at 
a radius of 2 Mpc, which is the edge of our detection limit. The SW merger significantly boosts 
the surface brightness, allowing us to detect X-ray emission out to as far as ~2.5 Mpc along this 
direction. The X-ray image also reveals two arc-shaped regions with excess surface brightness towards 
the east and west of the main cluster core. These regions appear over-pressured and they most likely 
originate from merger induced large scale supersonic gas motions, which are driving shocks into the 
ICM. The azimuthally averaged temperature profile, as well as the deprojected density and pressure 
profiles towards the E and NW, all show a sharp drop consistent with an outward propagating shock 
front located at the outermost edge of the giant radio halo observed at 352 MHz with the WSRT 
and which may be powering this radio emission. The entropy profiles along the relatively relaxed E 
and NW directions have a shape at large radii that is consistent with the average profiles of cool core 
clusters. This indicates that the accretion history (accretion rate, shock strength, and age) along the 
undisturbed directions of the non-cool core Coma Cluster are similar to those in evolved, well formed 
cool core clusters. Our data indicate a flat metal abundance profile at about 0.3 of the Solar value 
out to almost the virial radius of Coma, which would favor galactic winds over ram-pressure stripping 
as the dominant enrichment mechanism in the ICM. 
Subject headings: galaxy clusters 



1. INTRODUCTION 

All of the matter in the Universe, both luminous and 
dark, is distributed in a complex structure of sheets, fil- 
aments and voids which evolved following the gravita- 
tional collapse of small perturbations in the primordial 
density field. Clusters of galaxies are located at the nodes 
of this filamentary structure, also known as the cosmic 
web, and are constantly evolving and growing by accret- 
ing matter from the surrounding large-scale structure. 

The outskirts of galaxy clusters present an opportu- 
nity to study this process as it happens. The virial ra- 
diuf[^ approximately marks the border between regions 
of equilibration and in-fall. Merger activity from accret- 
ing material plays an important role in the dynamics 
of the outer regions. In order to understand large-scale 
structure formation in detail, it is therefore necessary to 
study the thermodynamic properties of the intra-cluster 
medium close to the virial radius in systems with a wide 
range of ages/masses and dynamical states. 

^ We adopt as our definition of the virial radius r200i within 
which the mean total density is 200 times the critical density of 
the universe at the redshift of the cluster 



Surface brightness profiles out t o the virial radius have 
been studied us ing ROSAT data ( [Vikhlinin et al.||1999 
Neumann 2005 ), which provided density estimates but no 
information about the ICM t emperature. Recen tly, using 



archival data from ROSAT, |Eckert et al. ( |2012[ ) find that 
beyond ~ rsoo, galaxy clusters deviate sigmhcantly from 
spherical symmetry, with only small differences between 
relaxed and disturbed systems. 

Robust measurements of the spectroscopic properties 
of the faint cluster outskirts have been made possible 
only recently, owing to the low and stable background 
of the Japanese-US Suzaku satellite. However, results 
to date have been obtained primaril y for m assive, re- 
laxed clusters such a s PKS 0745-191 (George ct al. 20091 
Walker et ar]2012a|), Abeh 2 204 (iReiprlch ct al. 2009L 
Abell 1795 ([Hautz et al.|2009| ,_AbeIljil3^( Hoshino^^ 



and 



2010|, the Ferseus Cluste r ([Simionescu et ai.||2011| )" 
XEell 2029 (jWalker et al.||2012c|. By comparison, fewer 



results are available tor less massive and/or less dynam- 
ically relaxed sys tems. For example, Abell 2142 (lAka^ 



matsu et al.|2011[ ) and Abell 3667 |Akamatsu et al.|20nf 
are so tar the only clear mergers for which the ther- 
modynamic profiles have been investigated out to r2oo- 



2 



Simionescu et al. 



Measurements of the properties of cl uster outskirts have 



also b een published for Abell 1689 (Kawaharada et al. 
2010| , t hought to be a line -of-sight merger, the Virgo 
Cluster (|Urban et al.||2011|), a young, dynamically un 



relaxed, low-mass system, the Hy dra A cluster, a low- 



mass relatively relaxed system ( jSato et al. 20121, and 
for the relaxed fossil gr oup /poor cluster KXJ 1159+5531 
( [Humphrey et al.||2012[ ). 



Here, we present results from a Suzaku Large Project 
targeted on the Coma Cluster, which is the nearest and 
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We assume 

17^=0.73, and Ho =70 km/s/Mpc. At th e redshift of the 
cluster, z=0.0231 ( Struble fc Rood|1999[ ), one arcminute 



corresponds to 28 kpc. 

2. OBSERVATIONS 
2.1. Suzaku 

A large mosaic of observations of the Coma Cluster 
was obtained as a Suzaku Large Project during AO-6, 
covering the E, NW and SW directions contiguously out 
to a radius of 2° (PI A. Simionescu). The mosaic consists 
of 24 pointings with clean exposure times between 5 and 
20 ks each, for a total observing time of 295 ks. In ad- 
dition, 6 more pointings were awarded during the same 
AO as a Guest Observer program (PI T. Sato), covering 
regions located towards the NE and W from the cluster 
center, with a total effective exposure of 175 ks. We com- 
bined these pointings with archival data obtained in AOs 
1-5, which totaled 580 ks of observing time. The details 
of all the pointings used in this work are summarized 
in Table [2] We label the pointings with letters indicat- 
ing the direction of the arm that the pointing belongs 
to, followed by numerical indices which increase with in- 
creasing radius, such that e.g. all pointings with index 4 
are offset by approximately the same distance from the 
cluster core. Pointings labeled with half-integer indices 
are marked such that for example 3.5 lies between and 
partly overlaps with pointings 3 and 4. The spatial cov- 
erage of our combined data set is shown in Figure [l]and 
comprises over 1.0 Ms of Suzaku time. 

The data were reduced using the tools available in the 
HEAsoft package (version 6.12) to create a set of cleaned 
event lists with hot or flickering pixels removed. All stan- 
dard recommended screening criteria were applied We 
only included observation periods with the geomagnetic 
cut-off rigidity COR > 6 GV. The X-ray surface bright- 
ness image, extracted in the 0.7-7 keV energy band and 
corrected for vignetting and instrumental background, is 
shown in the right panel of Figure IT] 

We used the latest calibration files which account for 
the modified non X-ray background of the XISl detec- 
tor following the increase in the charge injection level of 
2011 June 1; in addition, for the XISl spectral analy- 
sis, we excluded two columns adjacent on each side to 
the charge-injected columns (the standard is to exclude 
one column on either side). This is because the injected 
charge may leak into these additional columns and cause 

^ Arida. M^^ XIS Data Analysis, 

'http: / / heasarc.gsfc. nasa.gov/docs/suzaku/analysis/abc/node9. htrnl] 
(:^UiU) 



an increase in the instrumental background. We applied 
the latest contamination layer calibration from 2012 July 
19. 

2.2. ROSAT 

The Coma Cluster was also covered by a mosaic of 
four ROSAT Position Sensitive Proportional Counter 
(PSPC) observations extending out to radii of approx- 
imately 60' (with observation identifiers rp800006N00, 
rp800013N00, rp800005N00, and rp800009N00). These 
observations were performed between 1991 June 16-18, 
and their combined clean exposure time is 78 ks. The 
data were reduc ed using the Extende d Source Analysis 
Software (ESAS, [Snowden et ar]|1994[ ). Background and 
exposure corrected images in three energy bands (0.7-0.9, 
0.9-1.3, 1.3-2.0 keV respectively) were combined, hav- 
ing removed artifacts associated with the detector edges. 
The resulting image is shown in the top left panel of 
Figure [T] 

3. SPECTRAL MODELING 

We used the Suzaku mosaic to extract spectra from 
annuli centered on the cluster center. The projected 
and deprojected profiles of thermodynamic pr operties 
were obtained with the XSPEC (version 12.7, Arnaud 
19961 spectral fitting package, employing the modified 
C-statistic estimator. We used the projct model to de- 
project the data under the assumption of spherical sym- 
metry. Sets of spectra extracted from concentric annuli 
in the 0.7-7.0 keV band were modeled simultaneously. 
We modeled each shell as a single-temperature thermal 
plas ma in collisional ion ization equilibrium using the apec 
code ( Smith et al.|2001 ), with the temperature and spec- 
trum normalization as tree parameters. Unless otherwise 
noted, the metal abun dance was set to 0. 3 of th e So- 
lar value in the units of Grevesse & Sauval ( 1998 ) . The 
Galactic absorption column density was fixed to the av- 
erage value measured at the location of each respective 
Suzaku field from the Leiden- Argentine-Bonn radio HI 
survey ( Kalberla et al. 2005) ; the values used are given 
in TablelSl 

3.1. Background subtraction 

The cosmic X-ray background was modeled in regions 
free of cluster emission, using the four outermost point- 
ings in our Suzaku mosaic towards both the E and NW 
(which cover the area beyond 80 arcmin from the cluster 
core; yellow squares in Fig. [l]), as well as two ROSAT 
All-Sky Survey (RASS) spectra extracted from l°-radius 
circles centered 3° from the cluster core towards the E 
and NW. 

The background model consisted of a power-law emis- 
sion component that accounts for the unresolved popula- 
tion of point sources, one absorbed thermal plasma model 
for the Galactic Halo (GH) emission, and an unabsorbcd 
thermal plasma model for the Local Hot Bubble (LHB). 
Wc first fitted the Suzaku spectra in the hard band (2-7 
keV) and determined the power-law index, for which we 
obtain 1.50±0.07. This is slightly steeper tha n the widely 
used value of 1 .41 reported previously by e.g. |De Luca"^ 
Molendi| ( |2004| from XMM-Newton MOS data, although 
the ditterence is statistically not highly significant. 

Next, we fitted the Suzaku and RASS spectra in paral- 
lel, but using chi-square minimization for the RASS and 
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Fig. 1. — Left: ROSAT-PSPC image of the Coma Cluster in the 0.7—2 keV energy band; the locations of the Suzaku observations used 
for background estimation are shown as yellow squares; all other Suzaku mosaic pointings are shown as red squares. Right: X-ray surface 
brightness map of the Coma Cluster from the Suzaku mosaic in the 0.7—7 keV energy band. Point sources and other artifacts and structures 
excluded from the spectral analysis are shown in white. 



C-statistic for the Suzaku spectra. We used the 0.1-2 
keV band for the RASS spectra and 0.7-7 keV band for 
the Suzaku spectra. For the RASS data, we fi xed the 
power-law paramet ers to the values reported by |Kuntz 



Snowden 

T 



(2000) 
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1.46, Y = 8.88 X 10-^^Eo 
tons keV~^ cm~^ s~^ arcmin"^ at 1 keV) while for the 
Suzaku spectra, we fixed the power-law index to the value 
measured in the Suzaku hard band (F = 1.50) and left 
the normalization free. The temperatures and normal- 
izations of the GH and LHB components were left free 
in the fit, with the metallicities fixed to 0.3 and 1 solar, 
respectively. The results are summarized in Table [l] 

Using the 8 different Suzaku pointings, we are also able 
to estimate the spatial variability of the background com- 
ponents. We fitted each region separately with the model 
described above, allowing the normalizations of the GH 
and power-law components to vary, and estimated the 
standard deviation of these two quantities. We likewise 
report this in Table m 

An additional "hot foreground" component with 
kT~0.6-0.8 keV has been reported in the vicinity o f the 
Galactic pla ne, where it may be produced by stars (Ma- 
sui et al.||2009,). We do not expect this component to be 



important near the Coma Cluster, located at high Galac- 
tic latitude. Adding another component to the back- 
ground model described above, when fitting the RASS 
and Suzaku data in parallel, results in a detection of an 
~ 0.6 keV thermal emission that is formally statistically 
significant. However, this component is not required by 
the RASS data alone and likely represents a calibration 
mismatch between the two telescopes rather than a real 
thermal component. We therefore do not include such a 
component in our default background model. 



TABLE 1 

Cosmic X-ray Background parameters. Errors are given 
at the ac = 1 level. where applicable, the second error 
value represents the standard deviation between 
individual measurements in the different fields used. 





kT (keV) or T 


yt 


CXB 


1.50 


1.15 ±0.02 ±0.07 


GH 


0.19 ±0.01 


4.0± 0.7± 1.8 


LHB 


0.094 ± 0.005 


1.4± 0.1 



^spectrum normalization, J rienudV , in units of z)]'^ 
^ per 20'^7r sq arcmin area. 



cm 



Charge exchange between heavy ions in the solar wind 
and neutral atoms in the Earth's geocorona and in the 
solar magnetosphere produces emission lines of highly- 
ionized C, N, O, Ne, and Mg which can act as an addi- 
tional foreground to the cluster emission. The geocoro- 
nal solar wind charge exch ange (SWCX) flu x can vary on 
timescales of seconds (e.g. Fujimoto et aH2007) . Empir- 
ically, when the solar wind proton flux is below 4 x 10* 
cm~^ s~^, Suzaku spectra do not show strong SW CX 
signatures ([Fujimoto et al.|[2007 Yoshino et al. 2009). 

We calculated the solar-wmd jDroton flux at tne Earth 
using WIND SWE (Solar Wind Experimcnl]^ data and 
find that the proton flux reaches a level above 4 x 10* 
cm^^ s^^ during the observation dates of fields el, e2, 
e3, nw5, nw65, wl, w4, w5, w6, sw4. The maximum pro- 
ton flux during these observations is 8 x 10* cm~^ s~^, 
which means we are dealing, if at all, with mild flares 
(during strong flares, the proton flux can reach levels as 
much as an order of magnitude higher). The spectra 
of nw65 and w5 do not show anomalous features, or an 



3.2. Solar Wind Charge Exchange and Stray Light 



^ http:/ /web. mit.edu/afs/athena/org/s/space/www/wind. html 
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Fig. 2. — Surface brightness profiles in the 0.7-7 keV band ob- 
tained with Suzalcu towards the E (blue), NW (red), and SW (ma- 
genta). 

anomalous soft flux, compared to the seven Suzaku back- 
ground pointings during which the proton flux was low; 
moreover, the best-fit spectral parameters for the clus- 
ter emission agree well between the e3 and e35 pointings 
(the latter also shows a very low and stable proton flux) . 
Therefore we do not expect our results to be significantly 
influenced by SWCX flares. 

The effects of stray light contamination (light scattered 
from very bright sources outside the field of view, in par- 
ticular the cluster center) are mitigated by the choice of 
target - Coma is a non-cool core cluster, whose central 
surface brightness is an order of magnitude lower than 
in the Perseus Cluster, which has also been observed ex- 
tensively with Suzaku. We have found that stray light 
does not hinder robust measurements of the tempera- 
ture and density in the outskirts of the Perseus Cluster 
(Urban et al., in prep.), and thus do not expect signifi- 
cant effects in the case of the Coma Cluster where this 
effect should be an order of magnitude lower. This is 
supported by the fact that the best-fit spectral param- 
eters for the cluster emission agree well between the e3 
and e35 pointings - because of the mirror geometry, we 
expect that the two different roll angles of the e3 and 
e35 pointings would produce very different levels of stray 
light contamination, thus if stray light had an important 
effect, we should measure different parameters from the 
two observations. 

4. RESULTS 

4.1. Complex morphology of a merging cluster 

It is well known that the Coma Cluster is undergo- 
ing a merger, with a bright subcluster located towards 
the south-west of the main cluster center readily ap- 
parent in our ROSAT and Suzaku images. Based on 



Fig. 3. — ROSAT image divided by the az imuthally averaged 
profile . 352 MHz WSRT radio contours from [Brown fc Rudnick| 
i |2011| l are shown in black. 



XMM-Newton data, Neumann et al. (20011 argue that 



this southwestern substructure is experiencing its first 
infall into the Coma Cluster. The main cluster itself 
is a non-cool core system featuring two central galaxies 
with similar brightness, rather than a single dominant 
brightest cluster galaxy (BCG). Therefore, it is assumed 
that its formation history has included at least one other 
relatively recent major merger. This is confirmed by ob- 
servations in the radio band, which show that t he Coma 
Clust er hosts both a giant radio halo (Coma C, Willsonl 
1970 1 associated with merging activit y, and a peripheraT 



radio relic (1253-1-275, Jaffe & Rudnick 1979; Giovannini 
et al.|1985 ) in the direction of the infalling south-western 
subcluster. 

The ongoing merger with the south-western sub-cluster 
clearly affects the surface brightness distribution, which 
is significantly higher at radii between 25-80 arcmin 
along the merger axis compared to the other more relaxed 
directions observed with Suzaku (see Fig. [2]). Among the 
relatively undisturbed directions, the cluster is brighter 
towards the E out to about 30 arcmin, where the surface 
brightness drops sharply. An additional surface bright- 
ness bump followed by a drop is seen at ~40 arcmin. 
Along the NW, the surface brightness distribution also 
shows a bump followed by a relatively sharp drop at 
r ^ 40 arcmin. Interestingly, these surface brightness 
jumps coincide with the the outer edge of the bright ra- 
dio halo. 

To illustrate this, in Figure |3j we show the result of di- 
viding the ROSAT-PSPC image by the azimuthally aver- 
aged ROSAT surface brightness profile. Apart from the 
southwestern infalling subcluster, we note the presence 
of two arc-shaped regions with excess surface brightness 
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towards the east and west of the main cluster core. The 
presence of these structures was suggested by the inflex- 
ions in the surface brightness profiles described above, 
but Figure [3] shows a two-dimensional picture revealing 
the distribution and shape of these features. The outer 
edge of both of these X-ray bright regions coincides well 
with the outer edge of the giant radio halo observed at 
352 MHz with the Westerbork Synt hesis Radio Telescope 
(WSRT) ( [Brown fc Rudmck||20ll| ). In addition, an arc 
of enhanced surface brightness seems to track the entire 
northern edge of the giant radio halo, although at a much 
smaller contrast than the E and W features. 

4.2. Thermodynamic properties of the ICM 
4.2.1. Relatively Undisturbed Directions 

In Figure |4j we show the radial profiles of projected 
temperature and spectral normalization obtained from 
our spectral analysis for the relatively undisturbed di- 
rections covered by the Suzaku mosaic, namely the E, 
NE, NW, and W. We detect the Coma Cluster along 
these directions out to a radius of 70 arcmin (2 Mpc). 
Both the temperature and spectral normalization are in 
broad agreement across the arms and show consistent ra- 
dial trends among the different directions probed, with 
the temperature decreasing from about 8.5 kcV at the 
cluster center to about 2 keV in the outermost annulus. 
In the same figure, we also show the azimuthally aver- 
aged profiles obtained by fitting in parallel the spectra 
corresponding to the relatively undisturbed azimuths for 
a given radial bin, with a common temperature and spec- 
trum normalization. The temperature along the eastern 
arm is systematically lower than towards the NW and 
W from the center of the cluster out to 50 arcmin (1.4 
Mpc). Along the western direction, while the w5 point- 
ing located immediately outside 70 arcmin shows no sig- 
nificant emission above the background level, we detect 
a thermal emission with kT=0.5 ±0.2 keV and XSPEC 
normalization Y=0.51^Q^g (in the same units as Table 
[l|) in the w6 pointing which covers radii from approxi- 
mately 95-110 arcmin (2.7-3.1 Mpc). Since the emission 
does not extend contiguously out to this radius (due to 
the lack of detection at 70-95 arcmin), it is likely asso- 
ciated with an in-falling group and not with the cluster 
itself. 

We furthermore performed a fit to the spectra in each 
annulus where the metallicity was allowed to vary as a 
free parameter. To improve the statistical precision, we 
tied all the metallicities corresponding to a given radial 
bin to the same value (the temperatures and spectral 
normalizations were still allowed to vary with azimuth for 
a given annulus). The result is shown in Figure [s] Apart 
from the last annulus (60-70 arcmin) and the cluster 
core, the best-fit values agree well with the value of 0.3 
Solar which was used above. 

4.2.2. The South- Western Merger 

In Figure [6j we show the projected temperature and 
spectral normalization profiles along the direction of the 
infalling southwestern subcluster. It is immediately clear 
that the ongoing merger has a major impact on the sur- 
face brightness distribution along this direction, which 
is significantly boosted compared to the relatively undis- 
turbed azimuths over a wide range of radii, from about 




radius (arcmin) 



Fig. 4. — Projected radial profiles of the temperature and spec- 
trum normalization obtained from the Suzaku data along azimuths 
not aligned with the infalling southwestern subcluster. The aver- 
age profiles and their la confidence intervals are shown with solid 
and dashed black lines, respectively. The spectrum normalization, 

/ UsnudV, is given in units of 4^jQ^(-]^_|_^)j2 cm~^ per 20^7r sq 
arcmin area. 




20 40 

radius (arcmin) 



Fig. 5. — Projected radial profile of the metallicity obtained from 
the Suzaku data along azimuths not aligned with the infalling 
southwestern subcluster. The horizontal line marks the value 0.3 
solar. 
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Fig. 6. — Projected radial profiles of the temperature and spec- 
trum normalization obtained from tlie Suzaku data along the 
southwestern direction, compared to other azimuths not aligned 
with the infalling subcluster. 

25 arcmin (700 kpc) to the edge of the detection hmit, 
with the peak excess located about 1.3 Mpc from the 
main cluster core. As a result of this boost, we detect 
X-ray emission out to much larger radii than towards the 
other directions. The statistical error bars are also signif- 
icantly shrunk because a very deep archival observation 
with an exposure time of 176 ks covers the radial range 
of about 63-78 arcmin. 

The projected temperature profile shows a steady de- 
cline from about 35-55 arcmin, in agreement with the 
temperature decrease seen along the more relaxed W 
and NW directions. Beyond this, the ICM tempera- 
ture remains practically constant over one Mpc in ra- 
dius (from about 55-78 arcmin). In the last annulus 
(78-88 arcmin), the best-fit temperature and spectrum 
normalization drop dramatically, which may indicate of 
the boundary of a shock. 

4.2.3. Deprojection Analysis 

For the E and NW arms, where we have contiguous 
radial coverage, we also produced deprojected profiles 
of the electron density, entropy, and pressure. These 
of course rely on the assumption of spherical symmetry, 
which is likely to be only a very rough approximation for 
an unrelaxed system like Coma - however, the similarity 
of the projected profiles in Figure [4] argues that this is 
still a sensible assumption. The results are shown in Fig- 
ure (TI The errors in this case have been derived using a 
Markov-Chain Monte Carlo (M CMC) analysis. 

|Planck Collaboration et "aL ( |2012j ) estimate the scale 
radms of the Coma Cluster to be ~ (47±1) arcmin 
by using the existing XMM-Newton mosaic to iteratively 
calculate the Yx = Mgas x kT parameter of the cluster 
and using the Yx — M^qq scaling relation calibrated from 
hydrostatic mass estimates in a nearby cluster sample 



Fig. 7. — Deprojected radial profiles of the electron density, en- 
tropy, and pressure obtained from the Suzaku data along the more 
relaxed NW and E directions. 



(Arnaud et al 
radii m l^'igure 



2010). We have thus chosen to scale our 
by this value of i?5oo- 
irection, inside a radius of ^30 arcmin. 



Along the E 

the electron density is systematically higher and the en- 
tropy is lower than toward the NW. The deprojected 
profiles indicate another density and pressure discontinu- 
ity at 40 arcmin. Even though departure from spherical 
symmetry may bias our deprojected thermodynamic pro- 
files, the discontinuities at ~ 30 arcmin and ~ 40 arcmin 
are also seen in the projected surface brightness profiles 
(Fig. [2]) and therefore are most likely real. 

Fitting the electron density profile with a power-law 
shape, Uf, ~ r~", we obtain best-fit indices ue = 
1.42 ± 0.13 and unw^ = 1-39 ± 0.08 using the full ra- 
dial range (6-70 arcmin, or 0.13-1.5 r^oo)- These fits are 
overplotted in Figure [7| The best-fit slope steepens sig- 
nificantly outside 30 arcmin (840 kpc, 0.6 rsoo) towards 
the NW {aNW,out = 2.22 ± 0.44); towards the east, the 
slope at large radii remains the same as the best-fit over 
the entire radial range {aE,out = 1-40 ± 0.63), although 
this may be affected by the ringing in the deprojected 
profile. 

Also in Figure [71 we compare our measure ments of the 
ICM pressure wi t h the results presented by iPlanck Col-| 
laboration et al. (2012). We find a good agreement be- 
tween the buzaku measurements and the best fit model 
to the azimuthally averaged Pl anck data (Model C in 
Planck Collaboration et al.|2012 ). We caution the reader, 
however, that the data used to determine the Planck 
model includes the direction of the SW merger, whose 
thermodynamic properties we show to be markedly dif- 
ferent than the more relaxed NW and E directions for 
which we have calculated the deprojected pressure from 
the Suzaku data. The agreement between Planck and 
Suzaku therefore may be partly a coincidence, since we 
would expect the Planck pressure profile excluding the 
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SW direction to be lower than the Planck radial average 
over the entire cluster and thus potentially lower than the 
pressure measured from the X-rays for the same region 
of sky. 

5. DISCUSSION 

5.1. Entropy at large radii 

Standard large scale structure formation models show 
that matter is shock heated as it falls into clusters under 
the pull of gravity. Simple theoretical models of this 
process predict that the entropy K should behav e as a 
power-law with radius, K (x , w ith (3 ~ 1.1 — 1.2 (Tozzi 



fc Norman|2"00T||Voit et al.|20 05'). However, as more and 



more measurements ot the thermodynamic profiles of the 
ICM at large radii are being performed, primarily with 
the Suzaku satellite, it is becoming apparent that the 
entropy profiles flatten away from this expe cted behavior 
near th e virial radius. This has motivated I Walker et al.l 
( 2012b I to propose a new functional form to describe the 
observed entropy profiles in the ICM, namely: 

5/5(0.3r2oo) = A(r/r2oo)^-^e-('^/^'--'')' 

Assuming that r2oo = 1-5?'50Qj we fitted our observed 
entropy profiles in the Coma Cluster with this functional 
form and obtain A=4.71±0.46 and B=1.18±0.28. This is 
in agreement with the parameters determined by [Walker 



et al 
with 



(|2012b|), who found A=4.4t°:5 and B=1.0+^;^^, and 
the parameters found for the relaxed directions of 
the Perseus Cluster, A=4.75±0.16, B=1.06±0.09 (Urban 
et al., in prep.). We over-plot both the expected power- 
law shape with index 1.1 and the best- fit functional form 
d escribed above in the e ntropy panel of Figure [7] 

'Lapi et al. (2010) and Cavaliere et al. ( 2011| predict a 
flattening of the entropy profiles in the oiitskirts of cool- 
core clusters because of the weakening of the strength of 
the accretion shock as clusters get older. In the frame- 
work of this interpretation, the consistency of the entropy 
profiles of the Coma Cluster with those found in cool core 
clusters indicates that the accretion rate and the strength 
of the accretion shock along the undisturbed directions 
of the non-cool core Coma Cluster are similar to those 
in evolved, well formed cool core clusters. The fact that 
the thermodynamic properties of the outskirts of cool 
core and non-cool core clusters are similar indicates that 
even though they developed along different evolutionary 
paths, with non-cool core clusters like Coma experiencing 
more violent mergers, today their accretion rate from the 
surrounding large scale structure, along relatively undis- 
turbed directions, proceeds at a similar rate. 

Alternatively, several other explanations for the fiat- 
tening of the entro py profiles in the clust er outskirts have 
been proposed. Simionescu et al. (2011) suggest that an 
inhomogeneous distribution ot the X-ray emitting gas at 
large radii can bias high the measurements of the gas 
density obtained from the bremsstrahlung emission, and 
thus bias the entropy measurements low. If this is the 
case, the similarity between the behavior of the entropy 
profiles measured in the Coma Cluster and in cool core 
clusters would imply a similar degree of clumping in the 
ICM at large radii, independent of the clusters' dynam- 
ical state. Numerical simulations currently predict that 
the gas clumping factor in the outskirts of mergi ng sys 



2012 1 due to a larger azimuthal variation of the gas den- 
sity and larger number of infalling substructures. How- 
ever, this effect could be compensated for if dense gas 
clumps were also destroyed more easily in more dynam- 
ically active systems; numerical simulations may not be 
able to reproduce the destruction rate of these clumps 
reliably due to the relatively low spatial resolution that 
is currently achievable. 



terns should be significantly boosted (e.g. |Vazza et al. 



Hoshino et al. (2010) and Akamatsu et al. (20111 on 
the other hand suggest that the entropy fiattenmg may 
be due to electron-ion nonequilibrium at large radii. If 
accretion shocks and merger shocks behave like classi- 
cal shocks, then they heat the ions practically instanta- 
neously, whereas electrons equilibrate with ions on a coUi- 
sional timescale. Within about a gigayear after the pass- 
ing of a shock front, therefore, the electron temperature 
near r2oo, which is measured from the bremsstrahlung 
emission spectrum, may not be representative of the true 
thermal energy content of the plasma, and this underesti- 
mation of the gas temperature could cause the flattening 
of the entropy profiles. If this is the case, the similar 
shapes of the entropy profiles in cool core clusters and 
in the dynamically active Coma Cluster would suggest a 
similar age and strength of the last shock that the gas in 
the outskirts experienced. 

5.2. Large scale gas motions, shocks, and the origin of 
the radio halo 

The radial entropy distribution toward the E shows a 
clear dip at r~ 30 arcmin. Within the statistical un- 
certainties of our measurement, the temperature profile 
appears constant across this feature. Therefore, the dip 
in the entropy profile appears to be mainly due to the 
enhanced density, and the feature is over-pressured with 
respect to the gas at larger radii. The excess X-ray emis- 
sion and overall lower temperature towards the E are 
most likely due to the infall of a small group centered 
around th e two galaxies NGC 49 21 an d NGC 4911, as 
argued by Vikhlinin et al. ( 1997 ) and Neumann et al. 
(|2003| . These galaxies are located at the same position 
as the E arc in Fig. [3| Interestingly, the outer edge of the 
E surface brightness excess corresponds to the edge of the 
giant radio halo observed at 352 MHz with the WSRT 
at this azimuth. Since the eastern arc is over-pressured, 
the infalling group is probably still moving supersoni- 
cally and should drive shocks into the surrounding ICM, 
which likely explains the second surface brightness edge 
seen at a larger radius of 40 arcmin towards the E. 

The NW profile shows a surface brightness jump at 
40 arcmin, corresponding to the edge of the radio halo, 
which extends furthest in radius along this direction. At 
the same radius, the deprojected thermodynamic profiles 
towards both the E and NW show a sharp drop in density 
and pressure. 

Along the western direction, [Planck Collaboration! 
et al. (2012) found a marked local steepening of the ra- 
dial gradient of the SZ signal measured with Planck at a 
radius of about 30 arcmin, and interpret this feature as 
a shock with a Mach number of about two. At the same 
radius, we observe an X-ray surface brightness excess in 
the form of the western arc in Fig. |3l whose edge coin- 
cides with that of the 352 MHz radionalo. However, the 
X-ray temperature profile along this direction does not 
drop at 30 arcmin, but is consistent with the presence of 
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a discontinuity further out, beyond a radius of 40 arcmin. 

While we do not detect significant temperature jumps 
in the shapes of the temperature profiles along individual 
arms, the azimuthally averaged temperature profile does 
seem to drop sharply beyond a radius of about 40 arcmin. 

We are therefore presented with a complex picture of 
the Coma Cluster ICM, with multiple edges which all 
appear to be over-pressured. Towards the E and W, the 
inner features, or "arcs" , at radii of about 30 arcmin, 
which are over-pressured due to a higher gas density but 
not due to a temperature jump, trace well the edge of 
the 352 MHz radio emission along these directions. Be- 
yond these features, at a larger radius of 40 arcmin, we 
find a steep temperature gradient which may represent 
an outward propagating shock front, possibly surround- 
ing the cluster center at all azimuths. This feature is 
associated both with a temperature and a density jump, 
and its radius corresponds to the largest extent of the 
spherically asymmetric giant radio halo (which stretches 
out furthest towards the N). 

Merger-driven shocks and turbulence ma y lead to the 
ampl ification of the magnetic fields (e.g. D olag et al. 
2002 1 and to (re )-acceleration of electrons (e.g. hinsslin &: 
Biermann 19981, giving rise to radio emission, i'he ob- 
served spatial coincidence between the observed jumps in 
the pressure profiles and the outer edge of the giant radio 
halo lend further support to the picture that supersonic 
large-scale motions produced by major mergers lead to 
the origin of giant radio halos. 

Along the direction of the ongoing SW merger, the 
ICM temperature shows a marked decline starting at 40 
arcmin, the same radius where we found a drop in the az- 
imuthally averaged temperature profile along the relaxed 
directions. This steep gradient is found in the radial 
range of 40-55 arcmin, beyond which the temperature re- 
mains practically constant over one Mpc in radius (from 
about 55-78 arcmin) . This could argue that the infalling 
southwestern subcluster generates mixing in the ICM at 
large radii which causes the temperature to homogenize 
over this wide radial range. In the last annulus (78-88 
arcmin), the best- fit temperature and spectrum normal- 
ization drop dramatically, which is indicative of the pres- 
ence of a strong shock. The presence of such a shock at 
this location could explain the powering of the SW radio 
relic through particle re-acceleration. Hints for the pres- 
ence of this s hock have also been report ed from XMM- 
Newton data (Ogrean & Briiggen 2012). However, the 
systematic uncertainties on this measurement are large, 
as shown in the appendix. While the uncertainties re- 
lated to the unresolved background AGN power-law flux 
can be mitigated by using the hard-band to fit for the 
normalization of this component in the last annulus, the 
uncertainties related to the GH component remain very 
large, preventing us from drawing a definitive conclusion 
regarding the presence of this shock. 

5.3. The chemical enrichment of the ICM 

Our data point towards a flat metal abundance proflle 
at about 0.3 of the Solar value out to almost the virial 
radius of Coma, wh ich is similar to what was found in 
the Perseus Cluster ( Simionescu et al.|2011 Urban et al. , 
in prep.). In the last annulus, by contrast, the 2a upper 
limit on the metallicity is only 0.16 of the Solar value. 
This apparent drop in metallicity may result from multi- 



temperature structure in the ICM, which can cause the 
abundance measurement to be biased low especially at 
lower ICM temperatures, simila r to what was fou nd at 
large radii in the Virgo Cluster ( Urban et al.||20lT ). 

The two main mechanisms enriching the ICM with 
metals are ram- pressure stripping of galaxies moving 



through the ICM (|Gunri fc Gottll972||Nulsen et al.|1982 
Roediger & Briiggen " 2007[ 2008p and galactic winds 



driv en by the energ y of a large number of supernovae 
(e.g. Heckmaii||2003 ) . Because ram-pressure stripping is 
proportional to the iCM density and to the square of the 
relative velocity between the galaxy and the ICM, both 
of which increase with decreasing distance to the cluster 
center, it is expected to produce a metallicity profile with 
a negative gradient all the way to large radii. Urban et al. 
(in prep.) use the fact that the observed metallicity pro- 
files are flat to conclude that the dominant enrichment 
mechanism are galactic winds, which were particularly 
strong and dominated the enrichment of the intergalac- 
tic medium at z = 2 — 3. The ICM observed at large radii 
was thus most likely pre-enriched before it fell into the 
gravitational potential of clusters. Our measurements of 
the metallicity profile in the Coma Cluster support this 
interpretation. 

6. CONCLUSIONS 

We have presented results from a Suzaku Large Project 
targeted on the Coma Cluster, which is the nearest and 
X-ray brightest hot (^-^8 keV), dynamically active, non- 
cool core system, focusing on the thermodynamic prop- 
erties of the ICM on large scales. 

Along the relatively undisturbed directions covered by 
the Suzaku mosaic, namely the E, NE, NW, and W, we 
detect emission from the ICM of the Coma Cluster out to 
a radius of 70 arcmin (2 Mpc). Both the temperature and 
spectral normalization are in broad agreement and show 
consistent radial trends along those directions, with the 
temperature decreasing from about 8.5 keV at the cluster 
center to about 2 keV at the virial radius. The temper- 
ature along the eastern arm is systematically lower than 
towards the NW and W from the center of the cluster 
out to 50 arcmin (1.4 Mpc). Towards the SW, an ongo- 
ing merger has a major impact on the surface brightness 
distribution, which is significantly boosted compared to 
the relatively undisturbed azimuths over a wide range of 
radii, from about 25 arcmin (700 kpc) to the edge of the 
detection limit. As a result of this boost, we detect X- 
ray emission out to larger radii than towards the other 
directions (as far as ~2.5 Mpc). 

Apart from the southwestern infalling subcluster, the 
surface brightness profiles show multiple edges around 
radii of 30-40 arcmin. The X-ray image reveals two arc- 
shaped regions with excess surface brightness towards the 
east and west of the main cluster core. These regions ap- 
pear over-pressured and they most likely originate from 
merger induced large scale supersonic gas motions, which 
are driving shocks into the ICM. While we do not de- 
tect significant temperature jumps in the shapes of the 
eastern or western temperature profiles alone, the az- 
imuthally averaged temperature profile, as well as the 
deprojected density and pressure profiles towards the E 
and NW, all show a sharp drop consistent with an out- 
ward propagating shock front located at 40 arcmin, cor- 
responding to the outermost edge of the giant radio halo 
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observed a t 352 MHz with the WSRT ( [Brown fc Rud- 
nick||2011[ ). The shock front may be powering this radio 
emission. 

The entropy profiles of the Coma Cluster along the 
relatively relaxed E and NW directions have a shape at 
large radii that is consistent with the average profiles of 
cool core clusters. This indicates that the accretion rate 
(which is also possibly related to the degree of gas clump- 
ing) and the strength and age of the accretion shock along 
the undisturbed directions of the non-cool core Coma 
Cluster are similar to those in evolved, well formed cool 
core clusters. 

Our data points towards a flat metal abundance pro- 
file at about 0.3 of the Solar value out to almost the 
virial radius of Coma. Because ram-pressure stripping 
should produce a radial gradient in metallicity at large 
radii which is not observed, we conclude that galactic 
winds must provide the dominant mechanism through 



which metals are stripped from galaxies and enrich the 
ICM/IGM. 
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TABLE 2 

SUZAKU OBSERVATIONS USED IN THIS WORK. 



Target 


ObsID 


Obs. date 


RA 


Dec 


Exposure 
(ks) 


riff 

(10^° cm-": 


PI 

) 


El 


ac806030010 


2011-06-16 


13 00 54.91 


-1-27 54 06.8 


9 


0.78 


Simionescu 


E2 


ac806031010 


2011-06-17 


13 02 07.06 


-1-27 53 47.8 


8 


0.80 


Simionescu 


E3 


ac806032010 


2011-06-17 


13 03 18.67 


-1-27 53 45.6 


5 


0.80 


Simionescu 


E3.5 


ao806032020 


2011-12-22 


13 03 43.20 


-1-27 54 27.4 


9 


0.79 


Simionescu 


E4 


ao806033010 


2011-06-17 


13 04 31.01 


-1-27 53 40.9 


19 


0.78 


Simionescu 


E5 


ao806034010 


2011-06-18 


13 05 42.86 


-1-27 53 25.4 


13 


0.99 


Simionescu 


E5.5 


ao80603S010 


2011-06-18 


13 06 26.76 


-1-27 53 39.8 


11 


0.99 


Simionescu 


E6 


ao806036010 


2011-06-19 


13 07 05.30 


-1-27 53 12.1 


10 


0.99 


Simionescu 


E6.5 


ao806037010 


2011-06-19 


13 07 48.62 


-1-27 53 08.5 


12 


0.99 


Simionescu 


E7 


ao806038010 


2011-06-20 


13 08 27.43 


-1-27 53 06.0 


14 


0.99 


Simionescu 


NWl 


ao806039010 


2011-06-20 


12 59 14.26 


-1-28 10 07.0 


7 


0.88 


Simionescu 


NWS 


ao806040010 


2011-06-20 


12 58 13.70 


-1-28 40 10.9 


9 


0.85 


Simionescu 


NW3.5 


ao806040020 


2011-12-22 


12 58 06.94 


-1-28 44 19.0 


17 


0.85 


Simionescu 


NW4 


ao806041010 


2011-06-20 


12 57 43.20 


-1-28 54 43.2 


19 


0.88 


Simionescu 


NWS 


ae806042010 


2011-06-21 


12 57 12.89 


-1-29 08 43.4 


18 


0.95 


Simionescu 


NWS. 5 


ao806043010 


2011-06-22 


12 56 56.21 


-1-29 16 51.6 


14 


0.97 


Simionescu 


NW6 


ae806044010 


2011-06-22 


12 56 37.08 


-1-29 25 05.9 


10 


1.00 


Simionescu 


NW6.S 


ao80604S010 


2011-06-22 


12 56 18.58 


-1-29 33 08.3 


16 


1.03 


Simionescu 


NW7 


ao806046010 


2011-06-23 


12 56 00.84 


-1-29 41 27.2 


13 


1.03 


Simionescu 


SWl 


ao806047010 


2011-06-23 


12 57 44.45 


-1-27 43 35.8 


7 


0.88 


Simionescu 


SWS.S 


ae806048010 


2011-06-23 


12 54 34.75 


-1-27 05 10.7 


15 


0.78 


Simionescu 


SW6 


ao806049010 


2011-06-24 


12 54 06.34 


-1-26 59 42.7 


11 


0.77 


Simionescu 


SW6.S 


ao806050010 


2011-06-24 


12 53 38.35 


-1-26 52 23.9 


16 


0.75 


Simionescu 


SW7 


ao806051010 


2011-06-2S 


12 53 10.46 


-1-26 46 26.0 


13 


0.75 


Simionescu 


NE4 


ac806021010 


2011-12-06 


13 03 24.84 


+28 30 40.3 


20 


0.93 


Sato 


NES 


ao806022010 


2011-12-06 


13 04 40.61 


-1-28 42 33.8 


26 


0.95 


Sato 


NE6 


ao806020010 


2011-12-05 


13 05 56.02 


-1-28 54 13.3 


40 


1.01 


Sato 


W4 


ao80602S010 


2011-12-09 


12 55 53.76 


-1-28 30 38.9 


17 


0.85 


Sato 


W5 


ao806024010 


2011-12-08 


12 54 39.07 


-1-28 42 39.2 


27 


0.94 


Sato 


W6 


ao806023010 


2011-12-07 


12 53 22.66 


-1-28 54 13.3 


45 


1.00 


Sato 


CO 


ac801097010 


2006-05-31 


12 59 42.41 


-1-27 54 22.0 


164 


0.91 


Sarazin 


WO 


ao801044010 


2006-05-30 


12 58 46.54 


-1-27 56 47.8 


73 


0.90 


Matsushita 


Wl 


aG802084010 


2007-06-21 


12 57 22.27 


-1-28 08 25.1 


29 


0.87 


Bregman 


NW2 


ao802082010 


2007-06-19 


12 58 31.32 


-1-28 23 38.0 


48 


0.85 


Bregman 


SW3 


ao802047010 


2007-12-02 


12 57 01.39 


-1-27 34 17.0 


25 


0.84 


Inoue 


SW4 


ao802048010 


2007-12-04 


12 56 06.02 


-1-27 25 30.7 


30 


0.84 


Inoue 


SW4.S 


ao803051010 


2008-12-23 


12 55 27.10 


-1-27 17 17.9 


176 


0.81 


Sarazin 


S4.5 


ae805079010 


2010-06-01 


12 57 33.43 


+26 55 34.0 


85 


0.74 


Irwin 



SYSTEMATIC UNCERTAINTIES 

We estimate the systematic uncertainties by varying the CXB power-law normahzation and GH normahzation by 
the standard deviations among the 8 Suzaku background fields analyzed (see Table [T]) . The results are presented in 
Figure |8] Most quantities are robustly measured - the largest effects are on the temperature and normalization of the 
last annuli towards the E and SW. 
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Fig. 8. — Systematic uncertainties due to the CXB variance estimated in Table[l] 



